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Abstract 

In a model of TeV right-handed (RH) neutrino by Krauss, Nasri, and Trodden, the sub-eV scale 
neutrino masses are generated via a 3-loop diagram with the vanishing see-saw mass forbidden by a 
discrete symmetry, and the TeV mass RH neutrino is simultaneously a novel candidate for the cold 
dark matter. However, we show that with a single RH neutrino it is not possible to generate two 
mass-square differences as required by the oscillation data. We extend the model by introducing 
one more TeV RH neutrino and show that it is possible to satisfy the oscillation pattern within the 
modified model. After studying in detail the constraints coming from the dark matter, lepton flavor 
violation and the muon anomalous magnetic moment, and the neutrinoless double beta decay, we 
explore the parameter space and derive predictions of the model. Finally, we study the production 
and decay signatures of the TeV RH neutrinos at TeV e^e~ colliders. 

PACS numbers: 
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I. INTRODUCTION 



One of the most natural way to generate a small neutrino mass is via the see-saw mech- 
anism P]. There are very heavy right-handed neutrinos, which are gauge singlets of the 
standard model (SM), and so they could have a large majorana mass Mr. After electroweak 
symmetry breaking, a Dirac mass term M/j between the right-handed and the left-handed 
neutrinos can be developed. Therefore, after diagonalizing the neutrino mass matrix, a small 
majorana mass ~ mj^/Mn for the left-handed neutrino is obtained. This is a very natural 
mechanism, provided that ~ 10^^~^^ GeV. One drawback of this scheme is that these 
right-handed neutrinos are too heavy to be produced at any terrestrial experiments. There- 
fore, phenomenologically there are not many channels to test the mechanism. Although it 
could be possible to get some hints from the neutrino masses and mixing, it is rather difficult 
to reconstruct the parameters of the right-handed neutrinos using the low energy data [3]. 

Another natural way to generate a small neutrino mass is via higher loop processes, e.g.. 
Zee model j^, with some lepton number violating couplings. However, these lepton number 
violating couplings are also subject to experimental constraints, e.g., fi 67, r — 67. In 
the Zee model, there are also extra scalars whose masses are of electroweak scale, and so 
can be observed at colliders i^. 

On the other hand, recent cosmological observations have established the concordance 
cosmological model where the present energy density consists of about 73% of cosmological 
constant (dark energy), 23% (non-baryonic) cold dark matter, and just 4% of baryons. To 
clarify the identity of the dark matter remains a prime open problem in cosmology and 
particle physics. Although quite a number of promising candidates have been proposed and 
investigated in detail, other possibilities can never be neglected. 

Recently, Krauss, Narsi, and Trodden considered an extension to the SM, similar to 
the Zee model, with two additional charged scalar singlets and a TeV right-handed neutrino. 
They showed that with an additional discrete symmetry the Dirac mass term between the 
left-handed and right-handed neutrinos are forbidden and thus avoiding the see-saw mass. 
Furthermore, the neutrino mass can only be generated at three loop level, and sub-eV 
neutrino masses can be obtained with the masses of the charged scalars and the right- 
handed neutrino of order of TeV. Phenomenologically, this model is interesting because 
the TeV right-handed neutrino can be produced at colliders and could be a dark matter 
candidate. 

In this work, we explore in details the phenomenology of the TeV right-handed (RH) 
neutrinos. We shall extend the analysis to three families of left-handed neutrinos and explore 
the region of the parameters that can accommodate the present oscillation data. In the 
course of our study, we found that the model in Ref. |3] with a single RH neutrino cannot 
explain the oscillation data, because it only gives one mass-square difference. We extend the 
model by adding another TeV RH neutrino, which is slightly heavier than the first one. We 
demonstrate that it is possible to accommodate the oscillation pattern. We also obtain the 
relic density of the RH neutrino, and discuss the possibility of detecting them if they form a 
substantial fraction of the dark matter. We also study the lepton number violating processes 
and the muon anomalous magnetic moment, and the production at leptonic colliders. In 
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particular, the pair production of NiN2,N2N2 at e~^e~/fi'^^~ colliders gives rise to very 
interesting signatures. The so produced will decay into Ni plus a pair of charged leptons 
inside the detector. Thus, the signature would be either one or two pairs of charged leptons 
plus a large missing energy. 

The organization is as follows. We describe the model in the next section. In Sec. Ill, 
we explore all the phenomenology associated with the TeV RH neutrino. In Sec. IV, we 
discuss the signatures in collider experiments. Section V is devoted to the conclusion. 

II. REVIEW OF THE MODEL 

The model considered in Ref. has two extra charged scalar singlets Si, S2 and a right- 
handed neutrino Nr. A discrete Z2 symmetry is imposed on the particles, such that all SM 
particles and 5*1 are even under Z2 but 5*2, Nji are odd under Z2. Therefore, the Dirac mass 
term L(j)Nji is forbidden, where is the SM Higgs boson. The see-saw mass is avoided. 

In the present work, we extend the model a bit further by adding the second TeV right- 
handed neutrino, which also has the odd Z2 parity. The reason for that is because with only 
1 TeV RH neutrino, it is impossible to obtain two mass-square differences, as required by 
the oscillation data. However, with two TeV RH neutrinos it is possible to accommodate 
two mass-square differences with the corresponding large mixing angles. We will explicitly 
show this result in the next section. The most general form for the interaction Lagrangian 
isi 

C = ULlClT2L^St +giaNiStLR + 92aN2StLR + V{SuS2) + h.C. 

+Mn^N^CNi + Mn,N^CN2 (1) 

where a,P denote the family indices, C is the charge-conjugation operator, and 1^(5*1, 5*2) 
contains a term Xs{SiS2)'^. Note that fajj is antisymmetric under interchange of the family 
indices. Note that even with the presence of the first term in the Lagrangian it cannot give 
rise to the one-loop Zee diagrams for neutrino mass generation, because there is no mixing 
term between the Zee charged scalar Si and the standard model Higgs doublet that can 
generate the charged lepton mass. 

If the masses of Ni, N2, Si, S2 are arranged such that M^-^ < M^^ < Ms^ < Ms^, Ni 
would be stable if the Z2 parity is maintained. The A^i could be a dark matter candidate 
provided that its interaction is weak enough. Also, Ni, N2 must be pair produced or produced 
associated with S2 because of the Z2 parity. The N2 so produced would decay into A'^i and a 
pair of charged leptons. The decay time may be long enough to produce a displaced vertex 
in the central detector. The 5*2 if produced would also decay into Ni, N2 and a charged 
lepton. We will discuss the phenomenology in details in the next section. 



In principle, there are terms like NiN24> and MN1N2. The latter explicitly gives a mixing between the 
two RH neutrinos, while the former also gives the mixing after the Higgs field develops a VEV. However, 
the mixing term can be rotated away by redefining the A'^i and N2 fields. Effectively, the Lagrangian has 
the form given in Eq. (^). 
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III. PHENOMENOLOGY 



A. Neutrino masses and mixings 

The goal here is to find the parameter space of the model in Eq. (P) such that the neutrino 
mass matrix so obtained can accommodate the maximal mixing for the atmospheric neutrino, 
the large mixing angle for the solar neutrino, and the small mixing angle for 6*13 

Amatm ~ 2.7 X 10^^ eV^ , sin^ 20atm = 1-0 
Amsoi ^ 7.1 X 10~^ eV^ , tan^ ^.^i = 0.45 , 

sin^ 2^13 < 0.1. (2) 

The three loop Feynman diagram that contributes to the neutrino mass matrix has been 
given in Ref. The neutrino mass matrix {M^)a/3 is given by 

{M^)a/3 ~ (47^2)3 ^sfapmi^gpQamtJ^p (3) 

where denote the fiavor of the neutrino. Note that in the Zee model, the neutrino 
mass matrix entries are proportional to fap such that only off-diagonal matrix elements are 
nonzero. It is well known that the Zee model gives bi-maximal mixings, which have some 
difficulties with the large- mixing angle solution of the solar neutrino [6|. Here in Eq. (0) 
we do not have the second Higgs doublet to give a mixing between the SM Higgs doublet 
and Si , and therefore the one- loop Zee-type diagrams are not possible. However, the mass 
matrix in Eq. ^ allows for nonzero diagonal elements, which may allow the departure from 
the bi-maximal mixings. 

The mixing matrix between fiavor eigenstates and mass eigenstates is given as 



^ C13C12 ■5i2Ci3 Si3 \ 

-S12C23 - S23S13C12 C23C12 - S23S13S12 S23C13 , (4) 

V ■523-512 ~ ■513C23C12 ~'523Cl2 — ■5l3'5l2C23 C23C13 / 

where we have ignored the phases. The mass eigenvalues are given by 

U^MU = Mdiag = diag(mi, m2, m^). (5) 
The mass-square differences and mixing angles are related to oscillation data by 



(6) 





= Am2i 


2 2 
= 1712 ~ ""^1 




= Am^2 


2 2 

= r/ig — m2 


dsol 


= di2 




a 

"atm 


= 023 ■ 





From Eq. Q the neutrino mass matrix is rewritten as 

A,, 



(47r2)3M5, 



^ (M^)e ifmg)eifmg)f, {fmg)eifmg)r 

{fmg)e{fmg)^ {fmg)l {fmg)f,{fmg)r \ , (7) 
\ {fmg)eifmg)r {fmg)^{fmg)r {fmg)l 
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where {fmg)a = J2p fapfnipdp, the mass eigenvalues are given by 

mi = 1712 = 0, (8) 



ifmg)l + {fmg)l + {fmg)l . (9) 



(47r2)3M52 

This model obviously cannot explain the neutrino oscillation data because of the vanishing 

Hereafter we would like to discuss a possibility to improve this shortcoming. The reason 
that this model predicts two vanishing mass eigenvalues is the proportionality relation in the 
mass matrix Therefore it is necessary to break the proportionality relation. Although 
one way to improve the mass matrix might be to add small perturbations to the original 
mass matrix, we however found that this approach cannot resolve the difficulty. Instead, 
we consider a modification of the right-handed neutrino sector. As mentioned before, we 
employ two TeV RH neutrinos, the mass matrix ((7j) is replaced by 

1 1 

(M^)„;3 ~ Xs ifmgi)a{9imf)/3, (10) 

where / denotes the two RH neutrinos. 

If we assume {fmg2)^ <^ {fmgi)^, Eq. (fTIH) is rewritten as 

/ 1 + w t + cd 



w w"^ wt I , (11) 

\t + cd wt + (f 



w = {fmgi)^/{fmgi)e, 
t = {fmgi)r/{fmgi)e, 

c= {fmg2)e/ifmgi)e, (12) 
d = {fmg2)r/ifmgi)e, 
and has one zero and two non-zero eigenvalues ; 

where 

2A± = l + w^ + t^ + c^ + d^ 



±^(1 + w^ + t^ + c^ + rf2)2 _ 4(rf2 + ^2^2 + ^2^2 _ 2cdt + cH^), (14) 

and each of the mixing angles is given by 

w{X+ - c2 - d^) 



^23 



t(A+ - c2) + cd 



(15) 



A+ — (i^ — ted 

Sl3 = I =, (16) 

^(A+ - rf2 _ tcdf + (1 + tl^)w\X+ - c2 - d^Y 

1 dw 

C12 = , (17) 

Cl3 J{c^ + rf2)^y2 + - df 
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where we adopt the normal mass hierarchy. Indeed, we found that the correct mixing angles 
could not be realized if we assumed the inverted mass hierarchy here. Here ^23 — 1, ■5^3 ^ 1 
imply w ^ t, A+ ^ c^, (P and w"^ S> 1. This means ~ ^ 1, c^, (P. Definitely, from 



tcd\ 



2 



,i„^29.3._h__ _ <0.1, (18) 



we obtain > 20. Since Eq. (fTTj) is rewritten as 



c^w"^ + {cw — d) 



2 



d^w^ 

where C13 ~ 1 and w ^ t are used, we obtain 

C2 1 
rf2 4' 

by comparing with Eq. From the mass-square differences ; 

A^L f^-V f2c^w^ + d^w^ -2cdwV 



we find 



^2c^w^ 


+ rf^w^ 




4^4 




2 


^2wy 


~ 10-' 


4 fw^ 




3 \ 20 





(20) 



(21) 



Finally, Am^^^ ~ m| = is rewritten as 

where we used A+ ~ 2^2. in the last subsection of this section, we find some parameter 
space that leads to correct mixing angles and mass-square differences, after considering also 
the constraints from the dark matter relic density and lepton flavor violation. 



B. Neutrinoless double beta decay 

A novel feature of the majorana neutrino is the existence of neutrinoless double beta 
decay, which essentially requires a nonzero entry (M,^)ee of the neutrino mass matrix. The 
nonobservation of it has put an upper bound on the size of {M^)^^ ^ 1 eV 0]. 

In the model with two RH neutrinos, {M^)ee is estimated to be 

{M.)ee - - (4^2)3V^J (/^^l)e + (M^2)^] 
/l ± 3 /20\\ 

^ 3 X 10-3 eV. (24) 

by using Eqs. (j^^ and ()23|) . Thus, we find that this model is consistent with the current 
experimental bound. Such a small (Mj,)ee jnay still be within the reach of the GENIUS 
neutrinoless double beta decay experiment 
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C. Dark matter: density and detection 



The lightest RH neutrino is stable because of the assumed discrete symmetry. Here we 
consider the relic density of the lightest RH neutrino, and the relic density must be less 
than the critical density of the Universe. First of all, we verify that the second lightest RH 
neutrino is of no relevance here because of the short decay time. The heavier RH neutrino 
will decay into the lighter one and two right-handed charged leptons, N2 {(^iP 
denote flavors), and its decay width is given by 



N2 



Mn, , ,2 1 



:|^l/3^2a| X 



2(1 - /is)(/il - /is)(/il + /is + /il/is - 3/ig) log 



5127r3'"^^"-' ' 2/^2 
+ (1 - /ii)/is(2/ii - 5/is - + 6/i^) - 2^1 log/ii 



/is - /^i 



(25) 



where /ii = M%^/ Mfj^, Hs = Mg^/M^^. In the worst case when Mj^^ is very close to Matj, 
say, they are both of order 1 TeV but differ by 1 GeV only, and we set gi ~ 0.1. In this 
case, the decay width is then of order 10"^ — 10^s~^, i.e., the decay time is still many orders 
smaller than the age of the present Universe. Therefore, the presence of N2 will not affect 
the relic density of A^i. 

The relevant interactions for the annihilation is NiNi ^tn^pR through charged scalar 
5*^ exchange. The corresponding invariant matrix element is given by 



\M\' 



\giagii3\ 



{2qi ■ pi)2q2 ■ P2 ^ {2q2 ■ Pi)2qi ■ P2 



(t - Ml ^2 



'52^ 



MI2? 



2Ml2pi-p2 
{t-MDiu-Ml: 



(26) 



where qi and pi are four-momenta of the incoming A^^i particles and the outgoing leptons, 
respectively. Then, we obtain 



2q^2q^av 



(Ppi d?p2 
(27r)22p? (27r)22p^ 

■M^J2 



[2'nY\M\'8^'\q, + q2-pi-p2) 



(27) 



Ml 



Stt [Ml + 

8 (M| - M2,J2 + |(M| - Ml) + ^ s /£ 

3 {Mi + i-M^y iU 



8 - f - - 



(28) 



where mia is the lepton mass. We expanded |7V1|2 in powers of the 3-momenta of these 
particles and integrated over the scattering angle in the second line. Following ref. the 
thermal averaged annihilation rate is estimated to be 



(av) 



'Ml T 




M, 



T 



T 



4(27r) 



9lagif3\ 



Ml + M^^ 



32n {Ml + Mi;)' 



■AM 



T 



AMlK, (^s/t) (2g°2g° 
( T \ 



av] 



Ni 



M. 



M 



(29) 
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where J2 / denotes the summation over lepton flavors, and we have omitted the contributions 
from the S-wave annihilation terms, which are suppressed by the masses of the final state 
leptons. The relic mass density is given by 



1.1 X 10' 



9 2{M^JT) 



g^Mp{av) 



GeV"^ 



(30) 



where is the decoupling temperature which is determined as 



In 



0.152 



- In In 

2 



0.152 



(31) 



and is the total number of relativistic degrees of freedom in the thermal bath (lo| . 



By comparing with the recent data from WMAP 
QnAjh^ = 0.113 = 2.2 x 10^^ 



11 



we find 



Uo3GeVy y^Mp^oM^, 



(32) 



We can calculate ctq from Eqs. (jHT|) and (1221), ^ind we obtain 

^2 \ 1/2 / 



1.4 X 10" 



10^ 



,g*iTd 



1 + 0.07 In 



( \ ( 10^ 

Uo^GeVy WiTd) 



GeV" 



(33) 



if we ignore the second term in Eq. ()31|) . Indeed, we can confirm the validity of this assump- 
tion within about 10% error by using Eq. (jHH|) . Actually, Eq. (j^ is evaluated to be 



ln(2.5 X 10^^) - -lnln(2.5 x 10^^) 
31 - 5.1 = 26. 



(34) 



Our result of {av) is consistent with a previous estimation [12]. Equations and (jHH)) 
read 



f 



M 



Ni 



1.3 X lO^GeV 



1 + 2' 



1 + 2 



1 + 



(35) 



It is obvious that the RH neutrino must be as light as ~ 10^ GeV and at least one of gia 
should be of order of unity, such that the relic density is consistent with the dark matter 
measurement. ^ As the mass difference between Ms^ and A^^i becomes larger, the upper 
bound on M^i becomes smaller provided that we keep g ^ 1. 

The detection of the RH neutrinos as a dark matter candidate depends on its annihilation 
cross section and its scattering cross section with nucleons. Conventional search of dark 



^ Krauss et al. |5|] claimed that Mjv^ ~ 1 TeV and ^ 0.1 is consistent with the dark matter constraint, 
but in their rough estimation a numerical factor of {To/Mm)/^ ~ 200 is missing from the equation of 

{(Tv). 
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matter employs an elastic scattering signal of the dark matter with the nucleons. We do 
not expect that the Nn dark matter would be easily identified by this method, given its 
very mild interaction. In addition, because of the majorana nature the annihilation into a 
pair charged lepton at the present velocity (frei ~ 0) is also highly suppressed by the small 
lepton mass, even in the case of the tau lepton. However, one possibility was pointed out by 
Baltz and Bergstrom 2| that the annihilation NiNi — > ^+£"7 would not suffer from helicity 
suppression. The rate of this process is approximately a/vr times the annihilation rate at 
the freeze-out. As will be indicated later, the dominant mode would be fi~^fi~'~f. There is a 
slight chance to observe the excess in positron, but however the energy spectrum is softened 
because of the cascade from the muon decay. However, the chance of observing the photon 
spectrum is somewhat better [l^ . 



D. Lepton flavor changing processes and g — 2 

There are two sources of lepton flavor violation in Eq. (0). The flrst one is from the 
interaction fai3L^CiT2Lj3S^ . This one is similar to the Zee model. (However, the present 
model would not give rise to neutrino mass terms in one loop because of the absence of the 
Si — (j) mixing.) The flavor violating amplitude of via an intermediate up would be 

proportional to \fa/3f/3p\- The second source is from the term giaNiS2^aR in the Lagrangian 
(^. The flavor violating amplitude of — > £^3 via an intermediate Ni would be proportional 
to IgiaQipl- We apply these two sources to the radiative decays of £^37 and the muon 
anomalous magnetic moment. 

The new contribution to the muon anomalous magnetic moment can be expressed as 



Aa, 



967r2 



/je I 



Ml 



+ 



Ml 



F^iMyMi: 



(36) 



where F2{x) = (1 — 6x + Sx"^ + 2x^ — 6x^lnx)/6(l 
x^O, and ^2(0.25) ^ 0.125. We naively put F2{x) -- 
we obtain 



- x)"^. The function ^2(0;) 1/6 for 
1/6 for a simple estimate. Therefore, 



Attf, = 3 X 10 



-10 



2' 



'2 X 102 GeV\ 



M 



S2 



(37) 



which implies that /23, /21, S'l/i, (72/^ can be as large as 0(1) for 0(200 GeV) Si, S2 without 
contributing in a signiflcant level to Aa^. 

Among the radiative decays fi —>■ e'j is the most constrained experimentally, B{ix ^ 
67) < 1.2 X 10^^^ The contribution of the our model is 



B{ji^ 67) 



av 
3847r 



I f^r /' 



liT J re I 



'iQ\giegip!\' 
Ml 



FliMljMl) + 



36|gf2e^2Mp 

Ml 



Fi{MljMl) 



(38) 
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where v = 246 GeV. Again we take F2{x) = 1/6 and 0(200 GeV) mass for 8^,82 for a 
simple estimate. 



B{fi 67) = 1.4 X 10 



-5 



4 

Ms, ) ' V Ms, 



(1/^ / px lO^GeV y , ^ ,2 X 10^ GeV ^ 



2 /2 X 10^ GeV^ 



'52 



< 1.2 X 10"^^ , (39) 



which imphes that 

\ferfr^,\ < 1 X 10"^ \gle9l^,\ < 1 X 10"^ |^2e^2/.| < 1 X 10^^ (40) 

This is in contrast to a work by Dicus et al. In their model, the couplings gis are much 
larger than /j/s. 

E. An example of consistent model parameters 

Here we summarize the constraints from previous subsections, and illustrate some allowed 
parameter space. The prime constraints come from neutrino oscillations. The maximal 
mixing and the mass square difference required in the atmospheric neutrino, and the small 
6*13 read 



fri^m^gi^ ~ fi^rmrgir > feiMm^gif, + ferrrirgir ^ \j Y ^ 10'^ Gev ) ^^^^ 

where the terms freT^eQie and f^et^eQie have been omitted because these terms are sup- 
pressed by electron mass. The large mixing angle and the mass square difference required 
in the solar neutrino are given by 

fTemeg2e + /rM"^M5'2/. ^ 2(/e^m^5(2/, + feTfnrg2T) > /Me"^e5'2e + ff,Tmrg2T, (42) 
f fe^,m^g2^, + fermrg2r \ U^) 

On the other hand, the dark matter constraint requires at least one of the gig., gi^, gir 
to be of order of unity. While the muon anomalous magnetic moment does not impose 
any strong constraints, lepton flavor violating processes, especially i?(yU — > ej) gives the 
following strong constraints 

\f,rfre\ < 1 X 10-^ , (44) 
Ifi'lefi'l^l, \92eg2fi\ ~ 1 X 10~^ . (45) 

Now, let's look for an example of consistent parameters. From Eq. ()4H) . we obtain 
\m^gi^\ ~ Im^girl, in other words l^fi^l > l^fi^l, and 

/tm > fe^ + fre- (46) 
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Since either gi^ or gie must be of order of unity from the dark matter constraint, we take 
gi^ ~ 1. From Eqs. (jl^ and with g2T — 0, we obtain 

/tm - 2/e/., \m^g2f,\ > \meg2e\, (47) 

and 

gl ^ 8/3 X lO-'gl ^ 0.27(^7i^/l)l (48) 
Equations (j46|) and (j47|) can be rewritten as 

l»^ + 7^, (49) 

where we find that a mild cancellation between /g^ and is necessary. For instance, 
fre/frti = —1/3. The stroug cancellation corresponds to the small 613. However, a cancella- 
tion with too high accuracy would require a which is too big by Eq. Therefore, one 
can say that this model predicts a relatively large mixing in 613. Now we obtain an example 
set of parameters which makes this model workable and they are 



Ifi'lel ~ 


1 X lO"^ 


Ifi'lel 


-1, 


\9it\ 


~ 0.06, 


lfl'2e| ^ 


2 X lO"^ 




~ 0.5, 


\92t\ 


< 10^2 



/e^ ~ 1 X 10-^ U^2x 10-2, ^ -/,^. (50) 



IV. PRODUCTION AT e+e COLLIDERS 



The decay of N2 may have an interesting signature, a displaced vertex, in colliders. 
Depending on the parameters, N2 could be able to travel a typical distance, e.g. mm, in 
the detector without depositing any kinetic energy, and suddenly decay into Ni and two 
charged leptons. The signature is very striking. 

The NiNi, N2N2 and N1N2 pair can be directly produced at e^e^ colliders. The differ- 



ential cross section for e'^e 
da , , 



NiNi, 1 = 1,2, is given by 



dcos9 



[e e 



NiNj) 



Pi 



9fe 



2567r s 



Ml? 



it - MlY 



+ 



Ml,? 



[u - Mlf 



{t-Ml){u-Mi: 



Ml^ - s/2(l -/3/cos^), u 



where = ^1 - AM^Js, t 
total cross section is obtained by integrating over the angle 9: 

2{xj - - ^2 



Ml^ - s/2{l + /5/cos( 



(51) 
The 



NiNi 



9je 



+ Xs 



647rs -2xf + xj{6xs + 1) - 2xix,{3xs + 2)+ Xs{l + x,){l + 2xJ 



X 



Pi{-2xi + 2xs + 1) + 2{{xi - Xs? + Xs) log 



'2x1 -2xs + Pi -1 
2xi -2xs- Pi-l 



where xj 
given by 

da 



Mlfjs and Xg = Mg^/s. For N1N2 production the differential cross section is 



dcos9 



[e e 



NiN2 



\9le92e\ 

128n 



Pu 
s 



+ 



(n-M^J(n-M^J 
{u-Ml? 
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{t-Ml){u-Ml) 
and the integrated cross section is 



(52) 



cr(e+e" ^ N1N2 



\gieg2e\^ /?12 



1287r S 



/9l2s(-l + Xi + X2 - 2Xs){-l + /9l2 + Xi + X2 - 2Xs)(l + /3i2 - Xi - X2 + 2x5) 

X |/5i2s(— 1 + Xi + X2 — 2xs)(— 1 + I3I2 + 2xi — x\ + 2X2 — 6X1X2 — X2 
— 4xs + 8xiXs + 8x2X5 — 8Xg) 

+ s(2a/XiX2 + (Xi + X2)(xi + X2 - AXs - 1) + 2Xs(2Xs + 1)) 

/o2 f ^ , , o \2m /-I - /5i2 + + X2 - 2xA 1 

X (-1 + + - 2.,f) log j , (53) 



where j3i2 = \J {1 — Xi — — 4xiX2. The above cross section formulas are equally valid 
for collisions. Since the constraints from the last section restrict gi^. and g2e to be 

hopelessly small, we shall concentrate on using gi^ and g2^. 

The production cross sections for the N2N2 and N1N2 pairs are given in Figs, ^a) and 
(b), respectively, for y/s = 0.5, 1, 1.5 TeV and for M^Tj from 150 — 800 GeV, and we have set 
gifj_ = 1, g2fj. = 0.5 (see Eq. ljKnjl ). In the curve for N1N2, we set Mm^ = Matj — 50 GeV. We are 
particularly interested in the N1N2, N2N2 production, because of its interesting signature. 

As we have calculated the decay width of N2 in Eq. (j2SI), the N2 can decay into A^i plus 
two charged leptons, either promptly or after traveling a visible distance from the interaction 
point. It depends on the parameters involved, mainly the largest of gipg2a- As seen in Eq. 
(jKUjl the largest is |5'i^5'2/i| ~ 0.5, and so the decay of A^2 is prompt. Therefore, in the case 
of A''iA^2 production, the signature would be a pair of charged leptons plus missing energies, 
because the A^'i's would escape the detection. The charged lepton pair is likely to be on one 
side of the event. In case of N2N2 production, the signature would be two pairs of charged 
leptons with a large missing energy. Note that in the case of NiNi production, there are 
nothing in the final state that can be detected. From Fig. ^the production cross sections 
are of order 0(10 — 100) fb, which implies plenty of events with O(IOO) fb~^ luminosity. 



One may also consider S'2'S'2 pair production. The ^2 so produced will decay into S". 



^I'^aR ^2^aR^ where ia = e, /i, T. However, the constraints on the parameter space 
require the mass of substantially heavier than A^^i and A^2, and therefore the S2S2 pair 
production cross section is relatively much smaller. 



V. CONCLUSIONS 

In this paper, we have discussed a model that explains the small neutrino mass and dark 
matter in the universe at the same time. Such a model was proposed by Krauss et al. 
as a modification of Zee model. However, our study revealed that their original model is 
unfortunately not capable of explaining the neutrino oscillation pattern. 
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FIG. 1: Production cross sections for (a) N2N2 and (b) N1N2 pairs for ^/s = 0.5, 1.0, 1.5 TeV 
at collisions. We have set (71^ = 1, = 0.5, as suggested by Eq. (|50j) . = 500 GeV, and 
Mtvi = M7V2 - 50 GeV. 
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We have extended the model by introducing another right-handed neutrino. We succeed 
in showing that such an extension is possible to achieve the correct neutrino mixing pattern. 
A prediction of this model is the normal mass hierarchy. In addition, the undiscovered mixing 
angle ^13 is relatively large, because of the requirement of a mild cancellation between the 
parameters for a small ^13 and a sensible coupling of the charged scalar, A^. 

The relic density of the lightest right-handed neutrino has also been revisited. Under the 
constraint by WMAP we found that the mass of the right-handed neutrino cannot be as 
large as TeV but only of order 1 x 10^ GeV, after a careful treatment of the calculation. In 
addition, other constrains including the muon anomalous magnetic moment, radiative decay 
of muon, neutrinoless double beta decay have also been studied. With all the constraints 
we are still able to find a sensible region of parameter space. 

Finally, our improved model has an interesting signature at leptonic colliders via pair 
production of right-handed neutrinos, in particular N1N2 and N2N2. The N2 so produced 
will decay into A^i plus two charged leptons. Thus, the signature is either one or two pairs 
of charged leptons with a large missing energy. Hence, this model can be tested not only by 
neutrino experiments but also by collider experiments. 
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